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Abstract 

In this work, we investigate the diversity-muhiplexing tradeoff (DMT) of the multiple-antenna (MIMO) static 
half-duplex relay channel. A general expression is derived for the DMT upper bound, which can be achieved by a 
compress-and-forward protocol at the relay, under certain assumptions. The DMT expression is given as the solution 
of a minimization problem in general, and an explicit expression is found when the relay channel is symmetric 
in terms of number of antennas, i.e. the source and the destination have n antennas each, and the relay has m 
antennas. It is observed that the static half-duplex DMT matches the full-duplex DMT when the relay has a single 
^ antenna, and is strictly below the full-duplex DMT when the relay has multiple antennas. Besides, the derivation 

O of the upper bound involves a new asymptotic study of spherical integrals (that is, integrals with respect to the 

Haar measure on the unitary group U{n)), which is a topic of mathematical interest in itself. 

00 

P I. Introduction 

l-H 

c/3 Next generation wireless communication will take place over large networks with multiple-antenna (MIMO) 
lHj terminals. To be able to understand how to operate these complex networks optimally, we need a solid 
^ theoretical background. The diversity-multiplexing tradeoff (DMT) [1] provides such a theoretical back- 

j>. ground for multiple-antenna systems. It establishes the fundamental tradeoff between reliability and rate 

via diversity and multiplexing gains. Diversity gain is a measure of reliability, and shows how fast the 
^ error probability decays with increasing signal-to-noise ratio (SNR). Similarly, multiplexing gain is related 
to the transmission rate of the system, and shows how this rate increases with increasing SNR. 

^ Cooperation/relaying is going to be fully integrated in the standard operation of next generation wireless 

00 communications [2]. In wireless channels, nearby nodes can overhear source messages for free. Due 

O to this wireless broadcast advantage, relays can process the overheard information and forward it to the 

> destination terminal. The destination can then combine the direct signal from the source and the forwarded 

^ signals from the relays to improve the system performance [3], [4], [5]. 

To understand large, multiple-antenna cooperative networks, the theory has to account for practical system 
constraints such as power, bandwidth or delay limitations. Another important constraint in relay channels 
is the half-duplex constraint. As the transmit power overwhelms the received power, wireless devices 
cannot transmit and receive at the same time in the same band. In other words, they are not full-duplex; 
they have to operate in half-duplex mode. 

The half-duplex relay channel DMT has been already studied in [5], [6], [7], [8], [9] and [10]. In [5], 
orthogonal amplify-and-forward and decode-and-forward DMT are studied for the single-antenna relay 
channel. In [6], the authors study the single-antenna relay channel and investigate the DMT of the non- 
orthogonal amplify-and-forward and the dynamic decode-and-forward protocols. Prasad and Varanasi 
[7] investigate a multiple relay chaimel with single anteima source and relays and a multiple antenna 
destination and propose space-time coding strategies. Bletsas et al. [8] study relay selection as an alternate 

* Olivier LevSque is with the Faculty I&C, Ecole Polytechnique Fed6rale de Lausanne, 1015 Lausanne, Switzerland, olivier.leveque@epfl.ch 
** Christophe Vignat is with the Institut Caspar Monge, Universite de Mame-la-Vallee, 77454 Mame-la-Vallee, France, vignat@univ-nilv.fr 

* Melda Yiiksel is with the EEE Department, TOBB University of Economics and Technology, 06560 Ankara, Turkey, yuksel@etu.edu.tr 



2 



way to obtain DMT improvements in a multiple relay setting, and Pawar et al. describe in [9] a new 
"quantize-and-map" relaying scheme that achieves the half-duplex DMT. In all these works, the number 
of degrees of freedom in the source-destination channel is equal to one; i.e. either the source or the 
destination is allowed multiple antennas, but not both. In [10], the authors study the multiple- antenna 
relay channel DMT, thereby increasing the degrees of freedom in the direct link. They first state the static 
half-duplex DMT upper bound without explicitly computing it, and then show that the compress-and- 
forward protocol achieves this upper bound for any number of antennas at the source, the relay or the 
destination. 

The above mentioned papers enlighten many important issues about the DMT of the half-duplex relay 
channel. However, while compress-and-forward is known to be DMT achieving (under suitable assump- 
tions on the available channel state information), the exact static half-duplex DMT for arbitrary antenna 
configurations has not been computed yet. When nodes have multiple antennas, this computation is a 
mathematically challenging problem at first sight, as it requires the knowledge of the joint eigenvalue 
distribution of two correlated Wishart matrices, which has not yet been addressed in the literature. In this 
paper, we take a slight detour that allows us to express the static half-duplex DMT as the solution of a 
minimization problem. An explicit expression is moreover found for the DMT when the source and the 
destination have n antennas each and the relay has m antennas. This analysis brings us to the study of 



new asymptotic s of spherical integrals (see Section VI), which is a topic of mathematical interest in itself. 



Finally, a natural question regarding the half-duplex constraint on the relay is to understand whether this 
constraint imposes a limitation on the system performance, in comparison to a relay operating in full- 
duplex mode. The DMT computation leads us to the following rather unexpected answer: when the relay 
has two or more antennas, the half-duplex constraint imposes a strict limitation on the system performance, 
irrespective of the number of antennas at the source and the destination. The only case where the half- 
duplex constraint imposes no limitation (in the DMT sense) is therefore the case where the relay has 
a single antenna. This fact was already acknowledged in the case where source and destination have a 
single antenna each [10]; it is nevertheless surprising to observe that when the relay has multiple antennas, 
increasing the number of antennas at both the source and the destination does not allow to get rid of the 
half-duplex limitation. 

The paper is structured as follows. After the introduction of the system model in Section |II| we first review 



in Section III classical DMT results for the point-to-point MIMO channel and the full-duplex relay; we 
then derive the cut- set upper bound on the DMT in the half-duplex case and describe the DMT achieving 
compress-and-forward scheme. In Section |IV} we present our computation method for the DMT upper 
bound, that involves the study of the joint eigenvalue distribution of random matrices and leads us to 
writing the DMT as the solution of a minimization problem. In Section |V} we give the explicit expression 
of the DMT in the case where the source and the detination have the same number of antennas. Finally, 



we analyze in detail in Section |VI] the asymptotics of spherical integrals that are required for the DMT 
computation. 



II. System Model 

We study the MIMO static half-duplex relay channel where the source, the relay and and the destination 
have p, m and n antennas respectively, see Figure [1} The relay listens for a fraction t of the time, and 
transmits in the remaining (1 — t) fraction, t E [0, 1]. 

When the relay is listening, i.e. the system is in state qi, the received signals at the relay and the destination 
are 



^2,1 — HiXii + ^2,1 

^3,1 = GXii + Z^^i 
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Node 2, (m) 




Nodel,(p) Nodes, (n) 

Fig. 1. The half-duplex relay channel. The source, destination and relay have p, m and n antennas respectively. The relay listens for t 
fraction of the time, and transmits in the remaining 1 — t. 



When the relay is transmitting, i.e. the system is in state q2, the received signal at the destination is 

^3,2 = GXi^2 + -f^2-^2,2 + ^3,2 

The column vectors Xi^i and X2^u I = 1,2 are respectively of length n and m and denote the signals 
the source and the relay transmit in state qi. As the relay is half-duplex, X2,i = 0. Similarly F2,/ and 
Y^ i, I = 1,2 are respectively of length m and n and denote the received signals at the destination and 
at the relay. Note that 12,2 = 0. The channel gain matrices G, Hi, and H2 are of size n x p, m x p and 
n X m respectively and are assumed to have independent and identically distributed (i.i.d.) zero mean 
complex Gaussian entries with unit variance. Fading is supposed slow and frequency non- selective, so 
that it remains fixed for one frame length. The relay and destination noise vectors ^2,1 and Z^^i are of 
length m and n with i.i.d. complex Gaussian entries with zero mean and unit variance. 

Remark 1: If one considers the probably more realistic model where the variance of the entries of G, Hi 
and H2 are different for each matrix, the final result below concerning the diversity-multiplexing tradeoff 
does not change. 

The source and the relay have average short-term power constraints Pi and P2 for each codeword 
transmitted. As a constant scaling in transmitted power does not change the DMT [1], we will assume 
that Pi = P2 = P/2 in the rest of this paper. 

There is no transmitter channel state information at the source or at the relay. On the other hand, we 
assume that there are pilot signals to measure receiver channel gains, and thus the relay and the destination 
know their incoming fading levels. 

In this paper, we compute the static half-duplex DMT upper bound, when there is only receiver channel 
state information. On the other hand, the bound presented in the following is still valid, if transmitter 
channel state information is available at the source, and the relay and the destination know all the channel 
gains in the system, but the system operates at constant information rate. Under these assumptions, there 
exists a scheme that achieves the DMT upper bound we find. We describe this DMT optimal relaying 



protocol in Section III-B 



Although the relay is informed about its incoming channel gains, we assume that the relay does not use 
this information to determine the amount of time it listens or transmits. We therefore assume that the relay 
employs a static protocol for communication, in which the decision of transmitting or receiving is not 
based on the channel matrix realizations, but only on their respective distributions. Dynamic protocols, 
in which the relay chooses t based on its channel state observations, have the potential to achieve better 
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half-duplex DMT. However, dynamic protocols are more complex to study and the half-duplex DMT 
computation already constitutes a hard problem. In this paper, we thus focus on static protocols only. 

Without loss of generality, we also assume that the relay employs a fixed protocol. In fixed protocols, the 
relay does not carry information via breaking its transmission and reception intervals into smaller blocks 
and controlUng its state variable. If the relay protocol were not fixed but random, the state variable could 
be used to convey additional information. However, the increase is at most one bit and the DMT results 
remain the same [10]. 

As the source node does not have channel state information, an outage occurs if the mutual information at 
the destination is not large enough to support the fixed target communication rate the source chooses. In 
this paper, we study the minimum outage probability achieved by such a system at high SNR, following 
the approach of Zheng and Tse in [1]. In the next section, we overview the DMT for point-to-point MIMO 
channels, the DMT for the MIMO full-duplex relay channel, and the cut-set bound mutual information 
expressions for the MIMO static half-duplex relay channel. 



III. Preliminaries 

A. Diversity-Multiplexing Tradeoff for the Point-to-Point MIMO Channel and the Full-Duplex Relay 

The DMT was first defined for slow, frequency non- selective Rayleigh fading point-to-point MIMO 
channels [1]. Assume there are p antennas at the source and n antennas at the destination. Let R{P) 
denote the transmission rate of this n x p system and Pe(-P) denote the probability of error as a function 
of the transmit power (equal to the average received power as the channel gains have unit variance). We 
let the transmission rate vary logarithmically with the transmit power constraint. Then the multiplexing 
gain r and the corresponding diversity gain d are defined as 

lim f^^"^ ~ r and lim ^'^f _ 

P^oo log P P^oo log P 

The DMT curve dp^n{r) indicates the tradeoff between these two gains, and is the piecewise linear curve 
joining the points (A;, dp^n{k)) forO<A;<pAn = min{p, n}, where 

dp,n{k) ^ {p - k){n - k), k e {0, . . . ,p A n} 

Notice that rfp^„(r) = dn,p{r). For ease of notation, we will assume in the following that the transmission 
rate is chosen as R{P) — r log P, for all P. 

The full-duplex relay channel DMT was computed in [10]. The relay forms a virtual antenna array 
with either the transmitter or the receiver, so using the classical cut- set bound argument shows that the 
performance of the system is upper bounded by the minimum of the performance of a MIMO point-to-point 
channel with either {p-\-m) x n antennas or p x {m + n) antennas, i.e., 

dpuir) = min(rfp+„,„(r),rfp,„+„(r)) 

Moreover, the compress-and-forward protocol achieves this DMT upper bound, when the relay and the 
destination are informed about all the channel gains in the system [10]. In the following, we will use 
this performance as a benchmark and exhibit what loss is to be expected when operating the relay in 
half-duplex mode. 
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B. Achievable Scheme for the Half-Duplex Relay 

In this paper, our objective is to compute the static half-duplex MIMO relay channel DMT under the 
assumptions explained in Section [llj However, no relaying protocol is known to achieve this bound under 
these assumptions. Only under full channel state information assumption at the relay and the destination 
do we know how to achieve the static half-duplex DMT upper bound. In such systems, the only known 
static half-duplex DMT optimal protocol is compress-and-forward. 

In the compress-and-forward protocol, the relay first compresses its received signal 1^2,1, and then forwards 
the compressed signal through the relay-destination channel in the 1 — t fraction of the time it transmits. 
The compression at the relay is of Wyner- Ziv type [13], in the sense that the relay compresses its received 
signal taking into account that the destination has side information I3 1, available directly from the source. 
Although this protocol does not utilize the channel state information available at the source, the relay 
needs to know all the channel gains in the system to ensure that the compressed signal is received reliably 
at the destination. The soft information transmission of the compress-and-forward protocol is crucial to 
achieve optimal diversity and multiplexing gains. Protocols such as decode-and-forward [3], [11], which 
perform a hard decision about the source message, result in DMT suboptimal performance. 

As a side remark, we notice that the compress-and-forward operation is DMT optimal not only for the 
static half-duplex relay channel, but also for the dynamic one. However, for the dynamic half-duplex relay 
channel, the DMT upper bound is hard to compute. 

C. Cut-Set Upper Bound for the Half-Duplex Relay 

In [10], it is shown that cut-set bounds and the associated outage expressions can be used to find DMT 
upper bounds. Following the same approach, we provide next the cut-set bounds for the static half-duplex 
relay channel and the corresponding probability of outage expressions to find the best half-duplex relay 
channel DMT for static protocols. 

In the half-duplex, fixed and static relay channel, the mutual information expressions for the cut-sets 
around the source and the destination are respectively equal to 

hit) = tI{X,-Y^Y^\qi) + {l-mXv,Y:,\X^,q^) (1) 

lD{t) = t/(Xi;r3|gi) + (l-t)/(XiX2;F3|g2) (2) 

To find the best static DMT upper bound, we need to find the maximum of these mutual information 
expressions, which are obtained when Gaussian codebooks are used and the input covariance matrix is 
chosen optimally [11], [12]. 

We know that for any channel matrix H of size n x m and for any input covariance matrix Q of size 
m X m [1], 

sup log det (/„ + HQH*) < log det (/„ + PHH*) (3) 

Q>0,Tr{Q}<P 

where /„ denotes the identity matrix of size nxn, and * denotes conjugate transpose. Using this inequality, 
we can further upper bound ([T]) and ^ as Is{t) < I'sit) and /d(^) < Ini^)^ where 

^ + (l-t)logdet(J„ + PGG*) (4) 

I'oit) = thgdet (In + PGG*) + il-t)\ogdet (In + P[G,H2][G,H2Y) (5) 
Thus, the outage probability corresponding to a target rate r log P is lower bounded by 

Pout(rlogP) > min max{F{I's{t) < r log P), P(/^(t) < rlogP)} ^ Po„t,o(r logP) (6) 



I'sit) 



t log det Im+n + P 



G 




G 


Hi 




Hi 
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These expressions lead to the DMT upper bound 

d-HDir) < max mm{ds{r,t),dD{r,t)} = duDflir) (7) 

where Pout(r logP)=P-'^™(''\ ^{I'sit) < r logP)=p-'^s('^'*), Fil'^it) < r logP)=p-'^^('^'*) 
and Pout,o(rlogP)=P"'^™'°('') [] 

When the relay is half-duplex, it can only transmit during a fraction t E [0, 1] of the time (and therefore 
receive during the other fraction 1 — t). Because of the static protocol assumption, the fraction t is a fixed 
number, chosen according to the distribution of the channel coefficients only, and not to their realizations 
(but notice that t depends a priori on the target rate r). 



In Section IV below, we explain how to write the above diversity order dnDflir) as the solution of a 
minimization problem, following the methodology of [1]. In Section |V[ we show that in the particular 
case where the number of antennas at the source and the destination are equal (i.e., p = n), the symmetry 
of the above minimization problem allows us to write down (iHD,o(^) explicitly. 



IV. Computation Method 

In this section, our aim is to compute the diversity orders corresponding to the following outage proba- 
bilities 



Fil'sit) < rlogP) 



G 
Hi 



G 
Hi 



+ (1 - t) log det (/„ + PGG*) < r log P 



= P(^tlogdet (^/^+„ + P 
and 

P(/^(t) < rlogP) = P(t log det {In + PGG*) + (1 - t) log det (J„ + P[G,H2] [G,H2]*) < rlogP) 

for a given t G [0, 1]. As the dimensions p, m and n are arbitrary, it is sufficient to consider the second 
case; the first one will be deduced correspondingly. For ease of notation, let us also write H2 = H in the 
following, so 

Filoit) < rlogP) =F{t log det (J„ + PGG*) + (1 - t) log det (J„ + PGG* + PHH*) < r logP) 

where G, H are two independent matrices, each with i.i.d. A/c(0, 1) entries, G is n x p and if is n x m. 
Remember also that the diversity order is defined as 

, . ,^ log(P(/i,(t) <rlogP) 

dD{r,t) = - hm , r G 0,n (8) 

P-^oo logP 

Before entering into the detailed computation, let us make the following observation. Let Ai, . . . , A„ be the 
eigenvalues of GG* and z/i, . . . , i/„ be the eigenvalues of GG* + HH*. Since both these matrices are clas- 
sical Wishart matrices, the (separate) joint distributions of (Ai, . . . , A„) and (z/i, . . . , are well known. 
Nevertheless, in order to compute the above diversity order (following the methodology of Zheng and Tse 
in [1]), what is a priori needed is the joint distribution of all the eigenvalues (Ai, . . . , A„, i^i, . . . , z/„). This 
distribution is hard to obtain, because of the intricate correlation of the matrices GG* and GG* + HH*. 
A method is presented below that allows to compute the diversity order in the general case, but avoids 
the computation of the joint eigenvalue distribution. The latter question has not been addressed yet in the 
mathematical literature and remains an interesting open problem in its own right. 

'Note that f{P)=P~'^ means limp_oo log /(P)/ log P = c. Inequalities are defined similarly. 
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Our computation method goes as follows. We have 

t logdet(/„ + PGG*) + {l-t) logdet(J„ + PGG* + PHH*) 
= logdet(J„ + PGG*) + (1 - t) logdet(/„, + P//*(/„ + PGG*)-^H) 

Since G and H are independent and their respective distributions are unitarily invariant, this expression 
has the same distribution as 

logdet(/„ + PA) + (1 - t) logdet(/^ + PH*{In + PK)-^H) 

where A = diag(Ai, .... \n/\p, 0,...,0), Ai > ... > Xnhp are the non-zero eigenvalues of GG* and the 
number of zeros in the list is equal to (n — p)^ . Let now P = (/„ + PK)~^/'^H (n x m matrix) and 
observe that 

logdet(/„ + PH*{I^ + PA)-^H) = logdet(J™ + B*B) = logdet(4 + BB*) 
Conditioned on A (which is independent of H), the joint distribution of the entries of B is given by 

p{B\A) = det(7„ + PA)- exp(-IY(P*(7„ + PA)P)) 

The spectral decomposition of the n x n matrix BB* reads BB* = UMU*, where U is a n x n unitary 
matrix and M — diag(/ii, . . . , UnAm, 0, . . . , 0), where ni > ... > UnAm are the non-zero eigenvalues of 
BB* and the number of zeros in the list is equal to (n — m)+. 

As the diagonal matrices A and M might contain zeros, let us define A = diag(Ai, . . . , XnAp) as well as 
M — diag(//i, . . . , UnAm)- The Jacobian of the transformation B i— > (M, U) is given by 

nAm 

J{M, U) = (det M)\"\ A{Mf, where A(M) = Y[ {/ij - /ik) 

j<k 

Therefore, conditioned on A (or A), the joint distribution of (M, U) is given by 

p{M, U\A) = Cn,ra det(/„Ap + PA)™ exp(-Tr([/Mf/*(/„ + PA))) (det M)\"\ A{M)^ 

so U is independent of M and distributed according to the Haar measure on the set 14 (n) of unitary nxn 
matrices (that is, the columns of U form a set of n orthonormal vectors which are uniformly distributed 
on the sphere e C" : \z\ — 1}). 

This allows us to compute the conditional distribution p(M|A): 

p{M\A) = [ dUp{M,U\A) 

Ju{n) 

= Gn,m det{InAp + PA)" (det M)l™-"l A{Mf [ dU exp(-Tr(t/Mt/*(4 + PA))) 

Ju{n) 

= Cn,m det(4Ap + ^A)"" (det M)!"*""! A(M)2 exp(-Tr(M)) [ dU exp(-PTr(C/MC/*A)) 
Let finally 

nAp 

A(A) = n (^^- - ^'^^ 
j,k=i 

It is a well known fact that the distribution of A is the classical Wishart distribution: 

p(A) = Gn,p (det A)l"-fl A{Af exp(-Tr(A)) 
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SO we obtain 

p(A,M) = p(A)p(M|A) = C'„,^,pdet(J„;,p + PAr(detA)l"-^'l(detM)l"-"l 

xA(A)2 A(M)2 exp(-Tr(A + M)) / dU exp{-PTT{UMU*A)) (9) 



JU(n) 

Notice that the spherical integral on the far right may be rewritten as 

„ / / nAp nAm 

/ dU exp(-PTr([/Mf/*A)) = E^/ ( exp ( -P V V Xj Hk \ujk\ 

Mn) V V i=l k=l 

where E^/ denotes the expectation with respect to the Haar measure on U{n). 
Remember now that we are interested in computing 

Pout,o(rlogP) = P(t logdet(4 + PGG*) + (1 - t) logdet(J„ + PGG* + PHH*) < rlogP) 

p(A, M) dk dM 



' {A,M : log det(7„Ap+PA)+(l-t) log dct(7„Am+PA/)<r log P} 

Following Zheng and Tse, let us make the change of variables = P^"j and /i^ = P"^'=, where 
ai < . . . < ttnAp and Pi < . . . < PnAm- The behaviour in the limit P ^ oo of most terms in (|9]) is well 
known, except for the spherical integral. Let us define 



nAp nAm 



J(P) = I exp I - ^ 5^ P'-- 

3 = 1 k=l 



U.jk\ 



The asymptotic behaviour of X(P) is given by 

—oo, if 3j < n A p, k < n A m such that j + k = n + 1 and aj + < 1 



,„„Mi(p)) 

P^oo log P 



nAp nAm 

^ ^(1 — ctj — /Jfc)"'", Otherwise 
j=i k=i 

(10) 



We relegate the proof of this asymptotic equality to Section VI 



The diversity order corresponding to the cut around the destination ([8]) can now be computed via the 
standard method developed in [1]. The computation leads to the following optimization problem: 

nAp nAm 

dD{r,t) = min ^^(n + p — 2j + 1) + ^^(n + m — 2A; + 1) 

j=i k=i 

nAp nAp nAm 

-m - + J] - a . - (11) 

j=i j=i k=i 

where the minimization takes place over the set of variables q;„ap and PnAm > • • • > 

/?i > such that 

nAp nAm 

^(l-a,)+ + (l-t) 5^(l-/3.)+<r 
j=i k=i 

and aj + /Sa: > 1, for all j G {1, . . . , ri A p]. A; G {1, . . . , ri A m} such that j + k = n + 1. 
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Similarly, the diversity order corresponding to the cut around the source is given by 

pAra pAm 

ds{r,t) = min y^(j9 + n - 2j + 1) aj + y^(p + m - 2k + 1) Pk 

j=i k=i 

pAn pAn pAm 

j=i j=i k=i 

where the minimization takes place over the set of variables OpAn > • • • > > and PpAm > • • • > 
Pi>0 such that 

pAn pAm 

j=i k=i 

and aj + /3k > I, for all j G {1, . . . , p A n}, k E {1, . . . ,p A m} such that j + k = p + 1. 

In the next section, we solve the above minimizations problem in the particular case where n = p and 
thus find an explicit expression for the diversity order c?hd.o('^)- 



V. Explicit Expression for the Diversity-Multiplexing Tradeoff in the Case of Equal 
Number of Antennas at the Source and the Destination 

In the case where the number of antennas at the source and the destination are equal (i.e., p = n), the 
symmetry of the problem implies that the optimal value of Q a nd thus of ^ is reached when t = 1/2, 
for all values of r, and the minimization problems ( [TT] ) and p2[ ) boil down to the same expression: 

n nAm 

rfi3(r, 1/2) = 1/2) = min ^(2n - 2j + 1) a^- + ^(n + m - 2A; + 1) 

j=i k=i 

n n nAm 

-m 5^(1 - a,Y + 5^ 5^(1 - a, - (13) 

j=i j=i k=i 

where an > ■ ■ ■ > ai > and /?„Am > • . • > /3i > are such that 

n , nAm 

j=l k=l 

and aj + Pk > for all j G {1, . . . , n}, k E {1, . . . ,n A m} such that j + k = n + 1. So 

duBfiir) = dsir,l/2) = dD{r,l/2) 

The above minimization problem constitutes a convex programming problem. It therefore has a unique 
solution, which can be found via different classical algorithms. The solution given below has first been 
obtained through numerical simulations for a relatively large number of values of n and m (1 < n < 5 and 
1 < m < 12). From these results, the general form of the solution could then be deduced by extrapolation. 
We therefore acknowledge that our derivation is not completely analytical. Our confidence in the result is 
strong, however, and we check below that our proposed solution satisfies the constraints of the problem, 
so that it is at least an upper bound on the actual solution of the minimization problem. 

In order to describe the solution, let us define Iq to be the minimum of n and [^^^J (so Iq = n if 
m > 3/2 — 1; otherwise, Iq is an integer number between and n — 1). It turns out that this number 
delimitates three different regimes for the DMT curve: 
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d(r) 
30 , 




Fig. 2. Half-duplex DMT curve for n — 3 and m — 1 (left) / n — 3 and m = 3 (middle) / n — 3 and m = 7 (right). 



a) For < r < /o/2 (low multiplexing gain), the outage probability of the half-duplex relay is determined 
by the outage of the nxm relay link matrix H (the direct link matrix G being already off). In this regime, 
the comer points of the diversity curve are given by 

'^HD,o(V2) = + {m - l){n - I), where / G {0, . . . , /q} 

b) For lo/2 < r < n — Iq/2 (intermediate multiplexing gain; notice that this regime does not exist if 
Iq = n, that is, if m > 3n — 1), the outage probability is determined by a combination of direct link G 
and relay link H outages, and the corner points of the diversity curve are given by 

'^HD,o(^o/2 + 1) = ll + {n + m- l){n - /q - I), where / G {0, . . . , n - /q} 

c) For n — Iq/2 < r < n (high multiplexing gain), the outage probability is determined by the outage of 
the n X n direct link matrix G only, and the corner points of the diversity curve are given by 

duDfiin - 1/2) = P, where / G {0, . . . , /q} 

Finally, the diversity curve c?hd,o(^) is the (convex and) piecewise linear curve interpolating between all 
these corner points. 

Notice that, as already mentioned in the introduction, the only case where the half-duplex curve matches 
the full-duplex (m + n) x n curve is when /q = 0, i.e. m = 1, and n takes any integer value. For any 
m > 2, the half-duplex curve does not match (and is therefore strictly below) the full-duplex curve. This 
can be readily checked by noticing that the horizontal positions of the comer points do not match: in the 
half-duplex case (for m > 2), some comer points are located at non-integer multiplexing gains. 

The results are illustrated in a particular case on Figure |2j DMT curves are shown for the case where 
the source and the destination have n = 3 antennas each and the relay has m = 1, 3 and 7 antennas, 
respectively. As expected, the loss due to the half-duplex constraint increases as m gets larger. 

A precise description of the main outage events at the corner points of the diversity curve is given below, 
along with the computation of the corresponding diversity. The main outage events are described in terms 
of the number of links being active in the channel matrices G and H. For the direct channel matrix G, 
the number of active links is equal to the number of non-vanishing singular values a/Aj (or rank) of G, 
that is, the number of aj being equal to zero (since Xj = P~°j). For the relay channel matrix H, such a 
relation is less immediate, since the numbers Pj are related to the singular values yfjl] of the combined 
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channel matrix {In + PGG*)^^^'^H (through the equality fij = P^^^). From this relation, we can actually 
deduce that the number of active links in the relay channel matrix H is equal to the sum of the number 
of ttj and Pj being equal to zero. In more detail, this gives rise to the following interpretation. 

a) For r = 1/2, where / G {0, . . . , Iq}, the main outage event is the event that only / relay links are active, 
while all n direct links are off. More precisely, we have 

al = ... = a: = l and Pl = . . . = (3* = 0, (3^^ = . . . = p*^^^ = 1 

These satisfy the constraints and the terms in ([13]) become 



n mAn 

^{2n - 2j + 1) a* =n^, ^(n + m - 2j + 1) f3* = (m - l){n - /) 
i=i j=i 

n n,mAn 

j=l J,fe=l 

j + fc < ri 

so the corresponding diversity is given by 

c?HD,o(V2) = + (m - - I) 

b) For r = Iq/2 + 1, where / G {0,...,n — /q}, the main outage event is that only / direct links and Iq + / 
relay links are active. More precisely, we have 

"i = • • • = a? = 0, a^+i = . . . = a; = 1 and /5* = . . . = /5* = 0, p^^^ = ...= (3*^^^ = 1 

These satisfy the constraints and the terms in ([13]) become 



mAn 



^(2ri - 2j + 1) a* = {n - l)\ + m - 2i + 1) (3* = (m - h){n - k) 

j=i j=i 



n n,mAn 



m y [l-a*)+ =ml, ^ (1 - a* -/?*)+ = //q 



so the corresponding diversity is given by 



j=i j,fe=i 

j-{-k<.n 



c?HD,o(^o/2 + I) = (n + m ~ l){n - Iq ~ I) + I, 



c) For r = n — 1/2, where / G {0, . . . ,/o/2}, the main outage event is that only n — I direct links are 
active and all n relay links are active. More precisely, we have 

«* = ... = al_i = 0, = ... = «; = ! and pl = . . . = p* = 0, p*^^ = . . . = p*^^^ = 1 

These satisfy the constraints and the terms in ([13]) become 



n mAn 

J](2n - 2j + 1) a* = f, ^(n + m - 2j + 1) P* = {m-l){n- I) 
j=i i=i 

n n,mAn 

m ^(1 - = m{n - I), E " " = (n - 1)1 
so the corresponding diversity is given by 

rfHD,o(^ - V2) = 

Let us end this section by a comment on the more general case, where the number of antennas at the 
source is not equal to the number of antennas at the destination. In this situation, it is clear that due to the 
inherent asymmetry of the problem, the optimal proportion t where the relay is listening is not equal to 
1/2. In addition, this optimal proportion not only depends on t but also on the target rate r. This makes 
the computation of the diversity order highly cumbersome in this case. 



12 



VI. New asymptotics of spherical integrals 



The goal of this section is to prove ( [TO] ), namely that 

— oo, if 3j < n A p, k < n A m such that j + k = n + 1 and aj + Pk < ^ 



P^oo log P 

where 



nAp nAm 

^ ^(1 — ctj — /^fc)'*', Otherwise 
j=i k=i 



nAp nAm 

I{P) =Eu (exp[-J2Yl P^""''^' 

j=i k=i 



Ujkl 



Asymptotics of spherical integrals have been already studied in the mathematical literature (see for instance 
[14], [15]) and have moreover found applications to the study of the performance of MMSE receivers in 
CDMA systems, when users employ different powers (this can be inferred from the paper [16], although 
only the equal power case is considered in there). The asymptotic regime considered here is different 
from the one considered in the above references. In our case, the matrix sizes n, m are fixed and the 
parameter P tends to infinity, while in [14], [15], the matrix sizes n, m tend to infinity. Accordingly, our 
approach departs substantially from the above mentioned papers, since it highly relies on the fact that the 
dimensions n and m are both fixed, while P — > oo. 

Notice that in the case where the matrix dimensions are fixed, an alternate computation technique for 
spherical integrals was already developed in [17] in a slightly different context. In this work, the celebrated 
Harish-Chandra-Itzykson-Zuber formula was used in order to rewrite the integral as a determinant, and then 
the asymptotic analysis of the determinant was performed. This only works however when the dimensions 
p and m are greater than or equal to n. 

• Let us first show that if there exist 1 < j < n A p and 1 < k < n Am such that j + k = n + 1 and 
ttj + < 1, then 

(/ nAp nAm \ \ 

(~ ^ ^ pi-".-/^'^ \u,,\^^ j < exp(-P^) (14) 

where c = max{l — aj — Pk '■ l<j<nAp, l<k<nA'm, j + k = n + l}>0 (implying the first 
line of ([To])). Indeed, 

(/ nAp nAm \ \ / nAp nAm 

V j=i k=i J J \ ^"^ j=i k=i 

(nAp nAm \ 
-EE ^^■+'=."+1 ^'^"^^'''^ ^ exp(-P^) 
j=l k=l J 

where c > is defined as above. The equality (a) follows from the fact that if Ai > ... > XnAp and 
/"!>...> /inAm, thcu {ujk = Sj+k,n+i} minimizcs 

nAp nAm 

EE^j'^^i^j'^i^ 
j=i k=i 



over all the matrices in U[n). This shows (14). 

• Since we know that the integral decays super-polynomially with P if c > 0, let us now assume that 

ftj + /5fc > 1, VI < j < n A 1 < A; < n A m such that j + k = n + 1 (15) 
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Because of the ordering of the a's and /9's, this implies that aj + > 1 for all j, k such that j + k > n + 1. 
Let us also define 



riAp nAm 



j=i fc=i 

{j + k<n, aj+P^<l} 



It can first be checked that under the assumption ( |T5| ), 

J(P)=E(exp(-X(P))) 
Indeed, it holds under the same assumption that 

nAp nAm 

j=i k=i 

so 

E(exp(-X(P))) > J(P) > E(exp(-X(P)))e-''' 

• For each 1 <k < — 1) A m, let us now define jk = sup{l < j < {n — k) Ap : aj + (3k < !}• Notice 
that 

aj + < 1, VI < J < Jk 
because of the ordering of the a's, and also that 

jk < 3u HI <k 

because of the ordering of the /3's. Moreover, X{P) may be rewritten as 

(n-l)Am j^. 

x{p)= J2 



k=i j=i 

Next, we compute 



E(exp(-X(P))) = [ F{exp{-X{P)) >t)dt = [ P(X(P) < log(l/t)) rft 
Jo Jo 

From this, we deduce that 

1 

E(exp(-X(P))) > /"'P(X(P) < \og{l/t))dt > -P(X(P) < 1) 
Jo ^ 

and that for any < £ < 1, 

E(exp(-X(P))) < ^ + ^ P(X(P) < log(l A)) dt<e + P(X(P) < log(l/e)) (16) 
Remember that our aim is to show that E(exp(— X(P))) = P^'^, where 

nAp nAm (n-l)Am jf, 

j=l k=l k=l j=l 



It is therefore natural to choose e = P in (16), which leads to 

-P(X(P) < 1) < E(exp(-X(P))) < P"'^ + P(X(P) < rflogP) 

e 

If we then show that 

P(X(P) < 1) = P(X(P) < dlogP) = P-'^ (17) 
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we will have proved (fTOl). 



• At this stage, it is worth giving an intuition as to why ( 17 1 should hold, focusing on the case where 
m,p > n for simplicity. It can first be checked that 

P(X(P) < 1) = P < p-{i-.-/3.), VI < J < jk, l<k<n-l) 

Being part of an n x n unitary matrix U, the n"^ random variables IwjfcP are linked through n(n + l)/2 
constraints. Nevertheless, there are no more than n{n — l)/2 such random variables involved in the above 
expression, so the joint probability that they are all asymptotically small is not affected by these constraints, 
and we may as well consider them as independent in this limit. This leads us to the following estimate: 

n-l 3k n-1 jk 

p(x(p) < 1) = nn^(i"j''^i' - P'^'""'''^"^) = jj jjp-(i-°^-'3^-) 

k=l j=l k=l j=l 

which is the desired result (in the case where m,p > n). The rigorous proof follows. 

• Notice that since 

P(X(P) < 1) < P ( ^P^-"^-^*|Mjfcp < 1, VI < A; < (n- 1) Am I < P(X(P) < n) 



it holds that 



Jk 



P(X(P) < 1) = P I ^pl-":'~/5fc 



u 



■jk\ 



< 1, VI < A; < (n - 1) Am 



and similarly that 



3k 



P(X(P) < rflogP) = P I ^P^""^-^'=|Mjfcp < c/logP, VI < A; < (n- 1) Am 



Proving ( 17 1 therefore boils down to proving that 

/ 3k 



Ujk? < s, VI < A; < (n - 1) A m = P" 



(18) 



for s = either 1 or dlogP. This asymptotic equality may be recast into the following form: 



3k 



\Uik\ 



< 1, VI < A; < (n - 1) A m 



(n-l)Am jfc 

n n^^^ 

k=l 3=1 



.3=1 

where Sjk = either l/pi""j~/5fc or ((ilogP)/P^~"j~^'= is tending to zero as P ^ cxd, if j < j^. (notice that 
when present, the logP term in the s^^'s on the right-hand side does not affect the asymptotic polynomial 
decay in P). 



In order to estimate the above probability, let us expand it as 



3k 



ME 



\Ujk\ 



31 



< 1, VI < A; < - 1) Am 



Sji 



32 

< 1 I X P I ^ 

d=l 



Sj2 



< 1 



< 1 



3k 



ME 



\Ujk\' 
Sjk 



< 1 



< 1, Vl< / < A; 



X P 




< 1 



j,{n—l)Am 



E 



Sjl 



< 1, VI < / < (n- 1) Am 
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If we therefore prove that for all 1 < A; < (n — 1) A m, 

< 1 



■?fc I 12 
Uikl 



PIE'"' 



J2 ^ < 1, VI < / < A; 



■j=i -I 



jk 



(19) 



we will have proved (18 1. We will show the above inequality step by step, starting with the simplest case 
k = l. 

• The rest of the proof relies on the following fact (see for instance Lemma 2.2. in [18]): a matrix U 
which is Haar distributed on W(n) may be obtained from the Gram-Schmidt orthonormalization of an 
n X n matrix G with i.i.d. A/c(0, 1) entries, that is: 



Ml 


= Wi/ 




where 


Wi 


= 91 


U2 


= W2/ 


Ik2||, 


where 


W2 


= 92 


Uk 


= Wkl 


\\Wk\\, 


where 


Wk 


= 9k 



Uk,Wk,9k are the column vectors of the matrices U,W,G, respectively. A key observation is that condi- 
tioned on the column vectors ui, ...,Uk-i, the column vector Wk is a Gaussian vector with zero mean and 



covariance matrix /„ — Xlf=l^^^^^ 
• Step 1. First column (k=l). We need to compute 



•Jl I 12 
Will 



< 1 



0=1 



Sji 



ME 



I 12 



j=i J 



nE(^) i-..p< E kill 



where wu, . . . , Wni are i.i.d. ~ A/c(0, 1) random variables. For this, we need the following lemma, whose 
proof is relegated to the appendix. 

Lemma 1: Let A be a positive definite m x m matrix, Amin(^) denote its smallest eigenvalue and y ~ 
A/c(0, Im), z ~ A/c(0, /„) be independent Gaussian vectors. Then for all Aq > 0, there exist C2 > Ci > 
independent of the matrix A (but possibly dependent on the parameters m and n) such that if Aniin(^) > Aq, 
then 

JZl. < P(,M, < ) < J^. 

Remembering that the sji tend to zero as P — » 00, we can therefore apply the above lemma with 
A = diag f . . . , ^^^) to obtain 



' ( E (^) < E l-'nl' 

\i=i ^ Jl / j=ji+i 
which is the desired result for k = 1. 

• Step 2. Second column (k=2). We need to compute 



n 



Sji 



J2 



HE 



0=1 



\Uj2? 
Sj2 



< 1 



Jl I 12 



E 



< 1 



s,-i 



pipi v^<ik2ir 



Ml 



Jl I 12 
Will 



E 

i=i 



< 1 



Remember that conditioned on mi, W2 ~ A/c(0, /„ — miMi). Since rank(/„ — miMi) = n — 1, this covariance 
matrix may be rewritten as /„ — uiul = VV*, where V is an n x (n — 1) matrix; notice that the n x n 
matrix (mi|V^) is unitary, since Uiu\ + VV* = In, and that V*V = In-i- 
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The Gaussian vector W2 may also be rewritten as W2 = Vx2, where X2 ~ A/c(0,Jn-i). Noticing that 

—- x%y*Vx2 = \\x2\\'^, we obtain 



1^2 



J2 



0=1 



Sj2 



J2 



Ml 



ME 



0=1 



Sj2 



< 1^2! 



since the distribution of X2 does not depend on ui. Let now D = diag(si2^, . . . , Sj^2)^ ^ be the x — 1) 
matrix composed of the first j2 rows of V and tti be the vector ui restricted to its first j2 components. 



We see that VV* 



Ulul and that 



J2 



ME 



Sj2 



< l|3;2| 



f{x%V*DVX2 < \\X2\? 



Using the LQ decomposition of the matrix V (i.e. V = LQ, where L is a lower-triangular x (ri — 1) 
matrix and Q is a unitary — 1) x (n — 1) matrix) together with the fact that the distribution of X2 is 
unitarily invariant, we obtain 

p (x;v*DVx2 < \\x2f) = p {x;l*dlx2 < \\x2f) 



Writing now L = {L\0), where L is a square j2 x j2 matrix, and X2 = {y2,Z2), where 7/2 are the first j2 
components of X2, and Z2 are the last n — 1 — j2 ones (remember that n — 1 — j 2 > 1)^ we see that the 
above probability may be rewritten as 

P {xlL*DLx2 < llxaf ) = P (ylL* DLy2 < \\y2f + = P (vICl* DL - I,,)y2 < '4z^ 

Using again Lemma [T| we deduce that 

32 



¥[yl[L''DL-I,^)y2<zlz2 
Indeed, 



det{L*DL - Ij^) det{L*DL) det{D) det{LL*] 



det{LL*) = det{LL*) = det{VV*) G [1 - IImiH', 1] 

and the condition XljLi — ^ implies that l^jip < sji for all 1 < j < ji, in particular for all 
1 < j < j2, and therefore that 1 — ||mi|P is asymptotically close to 1. This concludes the proof for k = 2. 

• Step 3. k-th column. We repeat the procedure of Step 2 in this more general setting. We need to compute 



jk 



ME 



\Ujk[ 



< 1 



E 

\Wjk\ 
Sjk 



Sjl 



< 1, Vl< / < 



Ml, . . . ,Mfc-l 



E 



Sjl 



< 1, Vl< / < 



Remember that conditioned on wi, . . . , Uk-i, Wk ~ A/c(0, /„ — Yl^=i ^lU^). Since rank(/„ — Yl^=i = 
n — k + 1, this covariance matrix may be rewritten as /„ — Yl^=i = VV*, where ^ is an n x (ra — /c + 1) 
matrix; notice that the n x n matrix {ui\ . . . \uk^i\V) is unitary, since Yl'i=i "^i"^! + VV* = In, and that 

V*V = In-k+l- 



The Gaussian vector Wk may also be rewritten as Wk = Vxk, where Xk ~ A/c(0, /„„fc+i). Noticing that 



K'fcP = x*uV*Vxk = WxkW^, we obtain 



Jk 



ME 



Fife I 

Sjk 



< \\Wk\ 



Jk 



Ui, 



ME 

0=1 



\{Vxk) 

Sjk 



< \\Xk\ 
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since the distribution of Xk does not depend on ui, . . . ,Uk-i. Let now D = diag(s^^, . . . , V be 

the jk X (n — 1) matrix composed of the first jk rows of V and Ui, . . . , Uk-i be the vectors ui, . . . , Uk-i 
restricted to their first components. We see that VV* = Ij^ — Yl^Zl '^I'^l ^i^d that 




\{Vxk) 



" < \\xkf ) =f(^xlV*DVxk < \\xk\ 



Sjk 



Using the LQ decomposition of the matrix V (i.e. V = LQ, where L is a lower-triangular x — 1) 
matrix and Q is a unitary {n — 1) x (n — 1) matrix) together with the fact that the distribution of x^ is 
unitarily invariant, we obtain 

P (^xlV*DVxk < llxfcf ) = P [xlVDLxk < ||xfcf ) 

Writing now L = {L\0), where L is a square jk x matrix, and Xk = {yk, Zk), where yk are the first jk 
components of Xk, and Zk are the last n — k + 1 — jk ones (remember that n — A; + 1 — j/,. > 1), we see 
that the above probability may be rewritten as 

P {xlL*DLxk < \\xkf) = P {y*kL*DLyk < \\ykf + \\zkf) = P {y*kiL*DL - I.Jy^ < zlz^ 
Using again Lemma [T| we deduce that 



1.1 1 



P [L^DL-lAyk < zlzk) = , , = = ....n^ ....^^.^ =11 



det{L*DL - JjJ det{L*DL) det{D) dei{LL* 



Sj2 



Indeed, 

det(LL*) = det(LL*) = det{VV*) G 



fc-i 



1=1 



. I |2 

and the condition ^^Li < 1 for all 1 < / < A; implies that < Sji for all 1 < / < A;, 1 < j < j^, 
in particular for all 1 < j < jk, and therefore that 1 — J2^=i ll^dP asymptotically close to 1. This 



proves (19) and therefore concludes the proof of (10). 



VII. Conclusion and Perspectives 

In this work, we have expressed the half-duplex relay channel diversity-multiplexing tradeoff (DMT) as 
the solution of a minimization problem. In the symmetric case, i.e. when the source and the destination 
have n antennas each and the relay has m antennas, the solution to the minimization problem is given 
explicitly. Static relaying protocols are considered. Our results show that the static half-duplex DMT is 
equal to the full-duplex DMT only when the relay has a single antenna. In all other cases, performance 
losses are to be expected because of the half-duplex constraint. The case of the dynamic half-duplex relay 
is more complex and left open. 

We believe that the asymptotics of spherical integrals developed in this paper for the study of the half- 
duplex DMT are of interest in their own right, as such asymptotics were not considered before in the 



literature. Besides, as mentioned at the beginning of Section IV, the problem of determining the joint 



eigenvalue distribution of two correlated Wishart matrices remains open. 
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Appendix 

Proof of Lemma [7] First notice that 

/■oo 

¥{y*Ay < ll^f ) = / f>{y*Ay < t) /,i,|,2(t) dt 
Jo 

where /pj|2(t) = exp(— — 1)! and 

ny*Ay <t)= [ exp(-||yf ) dy = ^] [ exp{-CA-'0 

by the change of variable ^ = A^^'^y. This probability is further upper bounded by 

1 C f"^ 

^(y*^y ^ ^) ^ l{^ e : lief < t}| 



where Cm = \{^ e : < It is also lower bounded by 

F{y*Ay < t) > ] f exp (- ^ **^!' 1 d^ 

minimi / 

In total, we therefore obtain 

,2,/ Cm pt^'+^-i exp(-t) C™(m + n-l)! 1 



Tr^detAJo (n-1)! 7r'"(n-l)! det(A) 

and 



P(t/My<||^r) > " , / \„ """^ " ' dt 



TT™ det A Vl + Amm(^)/ Jo (n-1)! 
C„(m + n-l)! / K,UA) 1 



7r"^(n-l)! Vl + Vin(^)/ det(A) 
where we have used the change of variable s = ^ '^min(^^+^ j This proves the lemma, since by assumption, 

Amin(A) > Ao > 0. ■ 
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